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Summary 

We suggest several algebraic bounds for percolation on 
directed and undirected graphs: proliferation of strongly- 
connected clusters, proliferation of in- and out-clusters, 
and the transition associated with the number of giant 
components. 

Introduction 

Percolation on random graphs has been successfully used 
in network theory as a way to understand connectivity of 
real-life networks]!, 3, 13]. On a tree graph, mean field 
theory is useful to determine the location of a percolation 
transition [12]. However the existence of cycles reduce the 
applicability of this method on general graphs. We seek 
rigorous mean field bounds on the percolation transition 
which would be applicable to any graph. 

Spectra of associated matrices have been linked to per¬ 
colation and epidemic thresholds for a variety of graph 
and network models [10, 11, 4, 2, 9]. Recently, we 
constructed a lower bound on the percolation transi¬ 
tion for an infinite quasi-transitive graph a graph- 
theoretical analog of a translationally-invariant system[5]. 
The bound, 

Pc > 1 /p{H), (1) 

is defined by the inverse spectral radius of the Hashimoto 
matrix[6] H which is used to enumerate non-backtracking 
walks on Q q . The formal proof involved a sequential appli¬ 
cation of cycle unwrapping maps which results in a tree 
graph locally equivalent to the original graph. In the case 
of a degree-regular graph, Eq. (1) gives the well-known 
bound in terms of the maximum degree. In the case of 
random graphs with few short cycles, the same expression 
(1) gives a numerically exact result for the percolation 
transition [7]. 

In this work we analyze the applicability of the ex¬ 
pression (1) to more general graphs and in particular, fi¬ 
nite graphs. We note that the Hashimoto matrix H can 
be viewed as the adjacency matrix of the oriented line 
(di)graph (OLG) associated with the original graph [8]. 
Thus, Eq. (1) requires a mapping of the original undi¬ 


rected percolation problem to a directed one. Several tran¬ 
sitions are associated with percolation on a digraph T >, in 
particular, formation of a giant strongly connected compo¬ 
nent, and formation of a giant in- or out- component [11]. 
We show that on a general digraph, Eq. (1) is most di¬ 
rectly associated with another transition, the formation 
of a strongly connected component with a large number 
of distinct self-avoiding cycles (SACs), and argue that 
this property is related to uniqueness of the percolating 
cluster. 

Main results 

Consider an order-n strongly-connected digraph T> with 
vertex and edge sets V(£>), £(X>). Associated with V is 
an adjacency matrix A = A(T>), and Hashimoto matrix[6] 
H = H(V), with elements 

H-u,v — djj' (1 ^A), u = i y j■ v = j t f (2) 

where u , v are directed edges in £ (T>). The digraph T> may 
have some symmetric edges (length-two cycles), whereas 
the OLG has none. The spectral radius of H satisfies 
p(H) < p(A ); the equality is reached iff V has no symmet¬ 
ric edges. Also, for any induced g-norm, p(H) < ||-ffj| 9 . 

In site percolation on T> : each vertex is open with prob¬ 
ability p and closed with probability 1 — p; we consider 
a subgraph V of V induced by the open vertices. We 
are interested in the likelihood of forming in-, out-, or 
strongly connected components of a given size m on T>' 
[e.g., in the case of the out-component, there exists a ver¬ 
tex fo € V(D') such that m — 1 or more vertices can be 
reached by directed walks on V starting with fo]. In a 
strongly connected component, each site can be reached 
from any other by directed walks; it is automatically both 
an in- and out-cluster. We prove the following: 

Theorem 1 For any finite m > 0, the probability that a 
site v is the root of an out-component of size m or greater 
is bounded by mPm Ut \v) < (1 — p\\H ||i) _1 , p||ff||i < 1. 

This implies that out-cluster percolation threshold in 
large (di)graphs satisfies p£ out ^ > H-f/j]]" 1 . A similar state¬ 
ment can be written for in-cluster percolation, p c > 


1 


||.ff T ||j( 1 . This is a generalization of the maximum de¬ 
gree bound to arbitrary digraphs. The proof is based on 
constructing an upper bound for the average number of 
sites reachable from v via non-backtracking walks. This 
is found by upper-bounding the sum || Y^ m >o‘P rn H rn x v \\i, 
where x v has only one non-zero component. We con¬ 
structed a family of strongly-connected digraphs which 
saturate this bound, while at the same time p(H ) < ||i 2 j|i. 
Additionally, on these digraphs, the strongly-connected 
cluster percolation threshold is strictly higher than pc° ut ^ 
or pc m) ; they give a counterexample for the central con¬ 
jecture of Ref. [11]. 

We notice that OLG of a strongly connected digraph 
V is also strongly connected, provided that V remains 
strongly connected when any edge in any pair u = i —> j, 
u = j —» i of mutually opposite edges (if any) is removed. 
For such a digraph, an improved bound for out-cluster 
probability reads mPm Ut \v) < 71 , [1 — pp(-ff )] -1 (cf. Theo¬ 
rem 1 ), where 7 l is the principal ratio, 7 l = maxjj(£,/£.,•), 
for the left Perron-Frobenius vector of H, £,p(H) = fH. 
This inequality allows to extend the bound (1) from 
Ref. [5] to in- and out-cluster percolation on any strongly- 
connected quasi-transitive infinite digraph. 

The trace s -1 Tr H s counts non-backtracking directed 
cycles of length s on D. We constructed a corresponding 
bound for the total number of SACs on V, n E = \£(V)[. 

N < ^ s~ V Tr H s < n E |ln[l - pp{H)\ \ . (3) 

s> 0 

Notice that a strongly-connected cluster with m cycles 
supports m distinct SACs iff it has an articulation point 
separating any two cycles; such a cluster can be sepa¬ 
rated into two by removing a single site. A two-site- 
connected cluster (which can be cut into two pieces by 
removing two sites, but not one) will support at least 
m(m — l)/2 distinct SACs. A stable giant cluster which 
can not be easily separated into two large pieces will nec¬ 
essarily have large connectivity and support up to the 
maximum 0(2 m ) SACs. For a cluster occupying a finite 
fraction of directed edges in V , such a possibility is ex¬ 
cluded by Eq. (3). Therefore, we expect any large clusters 
formed at pp(H) < 1 to be unstable; i.e., fracture easily 
into two or more large cluster. Thus, we give 

Conjecture 2 For pp(H) < 1, if there is a percolating 
cluster, it is not unique with probability approaching one 
at large n. 


Conclusions 

We used algebraic techniques based on non-backtracking 
(Hashimoto) matrix to establish several new bounds for 
percolation on general directed (and undirected) graphs. 
Our results are formulated as bounds for the probability 
that a given site connects to a cluster of size m or greater, 
and are applicable for finite and infinite (di)graphs. We 
also discuss the stability of giant strongly-connected clus¬ 
ters, which is related to the uniqueness transition. 
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